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TECHNICAL NOTE 3067 

ROLLING EFFECTIVENESS AND AILERON REVERSAL OF 
RECTANGULAR WINGS AT SUPERSONIC SPEEDS 
By John M. Hedgepeth and Robert J. Kell 


SUMMARY 


linearized supersonic lifting-surface theory is used in conjunction 
with structural influence coefficients to formulate a method for ana- 
lyzing the aeroelastic behavior in roll at supersonic speeds of a rec- 
tangular wing mounted on a cylindrical body. Rolling effectiveness and 
aileron -reversal speed are computed by using a numerical solution which 
incorporates matrices. 

Results obtained for an example configuration by using this method 
are compared with the results obtained by using simplified methods of 
analysis. For the particular configuration considered, the variation of 
rolling effectiveness with Mach number is found for two constant-altitude 
flights. 


INTRODUCTION 


Tn the past, most aeroelastic calculations have been based on the 
use of beam theory for the structural analysis and strip theory for the 
aerodynamic analysis. The application of these simplified theories 
avoids complications which result from using more refined theories; in 
addition, the simplified theories are quite adequate, in most cases, 
when applied to wings of high aspect ratio. When applied to wings of 
low aspect ratio, however, these simple theories may become inadequate; 
if so, more refined structural and aerodynamic analyses are necessary. 

The purpose of this paper is to describe a method for predicting 
aeroelastic effects on the steady-state roll of rectangular wings at 
supersonic speeds in those cases for which beam and strip theory are 
inadequate but for which the aerodynamic effects of chordwise deforma- 
tion may be neglected. In this method the structural distortions caused 
by arbitrary loads are expressed in terms of structural influence func- 
tions. Tie aerodynamic loads caused by arbitrary angle -of -attack distri- 
butions are deter min ed by superposing basic aerodynamic loadings resulting 
from unit-step angle -of -attack distributions, which loadings are obtained 
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herein on the "basis of linearized aerodynamic theory. By means of this 
superposition procedure, first outlined "by Frick and Chubb in refer- 
ence 1, the application of three-dimensional lifting-surface theory is 
considerably simplified. 

The aircraft configuration considered herein consists of two flex- 
ible rectangular wings with trailing-edge ailerons of constant chord, 
diametrically mounted on an infinitely long, rigid, cylindrical fuselage. 

The analysis of the aeroelastic rolling behavior is separated into 
various parts. The analysis of the structural deformations is described, 
the aerodynamic loads are then found, and the two parts are combined. 

A numerical solution of the resulting equations is presented in matrix 
form. Tables of aerodynamic matrix elements usable for any rectangular 
plan form are included. A particular example is analyzed and the results 
are compared with those obtained by simplified methods. The structural 
analysis of the example configuration is included in appendix A and the 
details of the aerodynamic analysis are relegated to appendix B. 


SYMBOLS 


D 

E 

G 

GlCYjB) 


%(y^) 

L(y) 

M 


local flexural stiffness, Et^/l2(l - p^) 

Young’s modulus of elasticity 

shear modulus of elasticity, E/2(l + p) 

structural twist influence function which results from a 
unit concentrated load at the wing midchord 

structural twist influence function which results from a 
unit concentrated torque 

aerodynamic load per unit span, positive upward 

free -stream Mach number 


M(y) 


aerodynamic moment, per unit span, about the midchord, 
positive in the positive twist direction 


static pressure at altitude 


Pq standard static pressure at sea level 

aerodynamic moment, per unit span, about the elastic axis, 
positive in the positive twist direction 


Q(y) 
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V 

a 

b 

c 

c a 

e(y) 

l 

m 

P 

pb/2V 

q. 

t 

x,y,z 

P 

Pc z (y) 

P c m (y) 

P c q(y) 

6 

e(y) 

A 


free -stream velocity 

ratio between fuselage radius and exposed 'wing semi span 

total •wing span, 2(a l + l) 

single aileron span 

wing chord 

aileron chord 

distance measured forward from the midchord to the elastic 
axis, expressed as fraction of chord 

exposed wing semi span 

modified aspect-ratio parameter, pz/c 

rolling angular velocity (see fig. l) 

tangent of the wing-tip helix angle 

dynamic pressure 

thickness of wing cross section 

coordinate system (see fig. 1) 

cotangent of the Mach angle, ^M 2 - 1 


section lift coefficient. 


i»(y) 

qc/p 


section moment coefficient about the midchord. 


M(y) 

qc 2 /p 


section moment coefficient about the elastic axis, 

aileron deflection (see fig. l) 
angle of twist of wing (see fig. l) 


Q(y) 

qc 2 /p 


structural parameter, — \]2k(l - p) 
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|i Poisson’s ratio 

cr(x_,y) local angle of attack of wing 

0 ro lling effectiveness, (pb/2V) p/(pb/2V)^ 

’PCxjy) velocity potential 

$>(xi,yi) nondimensionalized velocity potential 

Subscripts: 

F flexible wing 

R rigid wing 

a airloads due to aileron deflection 

r airloads due to roll 

s airloads due to structural deformation 

rev aileron reversal 

p aerodynamic coefficients due to unit pb/2V 

p Q aerodynamic coefficients due to a unit rate of roll about 

the x-axis 

a aerodynamic coefficients due to a unit angle of attack of 

the entire wing 

8 aerodynamic coefficients due to a unit aileron deflection 

0 indicial aerodynamic coefficients due to a unit-step angle - 

of-attaek distribution 

I, II, ... VT regions on surface of wing 

1 nondimensional quantities used in appendix B 

Superscript: 

aerodynamic coefficients due to a unit-step angle -of -attack 
distribution on one wing only 


f 
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ANALYSIS 

Structural Deformations 


In this report the treatment of the structural deformations is 
based on the assumption that there is no chordwise bending of the wing. 
The effect of the structural distortion on the aerodynamic loads is 
then determined entirely by the angle of twist 9 (y) , the contribution 
of the spanwise bending being negligible. (See fig. 1.) Consequently, 
only the determination of 9(y) will be included in this analysis. 

If the section lift L(y) and the section moment about the mid- 
chord M(y) are known, the angle of twist can be obtained from 

e(y) = f Gl(y^) L (n) an + [ Gm(y^) m(ti) drj (1) 

'Jo <Jo 

In this equation, G-j- i (y, ri) and Gj^(y,T)) are influence functions which 

define the twist at y caused by the application of a unit concentrated 
load at the midchord and a unit concentrated torque, respectively, at 
the station q. For many structures these influence functions may be 
obtained analytically as is done in appendix A for a uniform flat-plate 
wing; other structures may be handled analytically by methods such as 
those described in references 2, 3, and 4, for instance. For Borne 
structures, it may be more convenient or even necessary to resort to 
the use of experimental influence coefficients. The analysis proceeds 
hereinafter on the assumption that the influence functions are known. 

In general, the two influence functions are needed in order to 
specify completely the twist of the wing. For many rectangular wings, 
however, sufficient accuracy can be obtained by expressing the twist 
solely in terms of the moment about some "elastic axis. " This elastic 
axis is herein defined as a line along which loads can be placed without 
producing significant twist anywhere. 

If an elastic axis does exist, it is no longer necessary to know 
the influence function associated with load; only Gj^(y,T]) need be 
determined. The twist, in this case, is given by 



6 


HACA TN 3067 


In this equation, Q(y) is the section torque about the elastic 
axis and can he expressed in terms of L(y) and M(y) by means of the 
simple moment-transfer relation 


Q(y) = M(y) - e(y) c L(y) (3) 

where e(y) is the distance measured forward from the midchord to the 
elastic axis, expressed as a fraction of the chord. 


Aerodynamic Loads 

The section lift L(y) and moment M(y) may be expressed in coef- 
ficient form as 


L(y) = -y pc z (y) 
M(y) = -- pc^y) 


00 


For convenience, the products of 3 = - 1 and the coefficients are 

considered herein rather than the coefficients themselves. 


It is assumed that the aircraft is undergoing a steady roll about 
the axis of the body and that this axis is in line with the direction 
of flight. Consequently, the resulting loads are due solely to struc- 
tural twist of the wing, the rolling velocity itself, and the deflection 
of the ailerons. Since linear aerodynamic theory is to be used, the 
coefficients resulting from this steady rolling maneuver can be expressed 
as 


P c z = (f*z ) 8 + (P c z) r + (e c z) a 

pc m = (P'mOg + (P'tiOj. + ( pc m) a 


(5) 


where s, r, and a refer, respectively, to structural deformation, 
roll, and aileron deflection. 

Tn the determination of the coefficients on the right-hand side of 
equations ( 5 ) the wing and body must be treated as a unit. In view of 
the fact that in this problem the pressures resulting from the presence 
of the body are important only in the neighborhood of the wing root and 
therefore contribute only slightly to the structural distortions and 
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rolling moment, a rather simple idealization is made regarding the 
bodyj that is, the hody is replaced hy a rigid plate joining the two 
wing roots together. The width of this plate is the same as the hody 
diameter. (See fig. 2(a).) While it Is not argued that this ideali- 
zation is the most nearly correct one, the resulting aeroelastic model 
does have the advantage of sinqplicity and allows aerodynamic interaction 
between the two wings. The presence of the cylindrical body is further 
taken into account by neglecting the effects upon the roll in g moment of 
the pressures acting on the rigid plate, since the pressures acting on 
the actual cylindrical body would produce no rolling moment. 

Although the details of the calculations of the aerodynamic loads 
have been relegated to appendix B, a short discussion of the loads due 
to each of the three causes - structural deformation, roll, and aileron 
deflection - is included in the subsequent portions of this section. 

Loads due to structural deformation . - As a consequence of the pre- 
viously assumed linearity of chordwise deformation, the local angle of 
attack of the wing is completely defined by the twist of the wing. Since 
the twist is not defined explicitly, the deformations being dependent on 
loads which are, in turn, dependent on the twist itself, it is necessary 
to be able to perform the rather difficult task of calculating the aero- 
dynamic loads caused by an arbitrary angle -of -attack distribution. This 
task is considerably simplified for the rolling problem by superposing 
loads caused by an antisymmetrical unit-step angle -of -attack distribution 
obtained by imposing a positive unit angle of attack outboard of any 
spanwise station T] on the right wing and a negative unit angle of 
attack over the corresponding portion of the left wing as shown in fig- 
ure 2(b). These basic load distributions, which have the nature of aero- 
dynamic influence functions, are hereinafter called "indicial" loads for 
the sake of brevity. Superposition of the indicial section coefficients 
of lift and twisting moment, designated 3cj 0 (y,q) and pc mo (y,T]) , 

respectively, which result from such basic angle -of -attack distributions, 
yields the section coefficients [^ c z(y)] B [P c m(y^ s due *’° s ' truc “ 

tural deformation. The required superposition integrals are given by 




pCjnCy) 

L- — *S 



cL0 

pc i0 (y ’’ l) Sf 4,1 




4,1 


(6) 


where the twist at the root is assumed to be zero. 
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The indicial section coefficients are dependent not only on the 
spanwlse coordinate y and the position of the step tj but also on 
two additional parameters. These parameters are the modified aspect- 
ratio parameter j3z/c, which usually appears in theoretical supersonic 
aerodynamic calculations, and the nondimens ional body radius a. The 
dependency on the body radius is undesirable since it restricts the 
application of these indicial section coefficients to a particular value 
of a. Fortunately the dependency can be elim i nated by makin g use of 
the fact that the loads due to the anti symmetrical unit-step angle-of- 
attack distribution can be separated into two parts: the first is the 

loads due to a unit-step angle of attack on the right wing only (see 
fig. 2(c)); the second is the loads due to a negative unit-step angle 
of attack on the left wing only. It can be seen that the second part 
is merely the negative mirror image of the first and, consequently, 
only the first case need be considered in detail. The total indicial 
section coefficients can be written in terms of these partial indicial 
section coefficients, designated pc^’ and fic^ ' , as follows: 


3c lQ (y,T]) = pc lQ '(y,ii) 
pCmgCy,^) = pcmg'Cy,!!) 


P c Z fl ' (-2aZ-y,Tj) 


3 c m 0 '(- 2al -y>n) 




(7) 


The first term on the right-hand side in equations (7) gives the contri- 
bution of the right-hand step angle of attack; the second term gives the 
contribution of the negative left-hand step. 


The advantage of the foregoing separation is that pc^' and pc m g ' 

are independent of the body radius and are functions of only the modi- 
fied aspect-ratio parameter, provided that this parameter is greater 

than — that is, when the Mach number is great enough so that 

there is no point on either wing that is influenced by both wing tips 

BZ 1 

simultaneously. Thus, if the restriction — > - — is imposed, a 

C JL "f* 

single parameter remains - pz/c itself - and it is feasible to compute 
tables applicable for any body radius. Numerical values of 3 c 2g' 


and pcjng* 


for - ~ % ^ % 1 and 

p z z 


0 -g 1 

L 


in steps of 0.1 are pre- 


sented in table I for several values of pz/c. These numerical values 
were obtained from expressions derived in appendix B. 
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Loads due to roll . - The section loading coefficients resulting 
from roll can be written in the form 


(e c 0 r = § 

( Pc m) r = 2V 


( 8 ) 


The rolling derivatives 


are given in appendix B. 


Pc lp and pc^ 

Here again the derivatives are dependent on both 3 l/c and a. Inspec- 
tion of the resulting expressions (eqs. (B15) and (Bl4) ) , however, 
reveals that the dependence on a is quite simplej the expressions 
for P c jp P c mp can eac ^ "b e separated into two parts - one with 

the coefficient a ' ’ ' ’ “ J ” ^ 


the other with 


- both of which parts 


1 + &' 1 + a 

are independent of a. It is possible to show that in each case the 
first part is merely the derivative which results from a unit angle of 
attack of the entire wing and the second part is the rolling derivative 
which would result if the wing were rolling about the x-axis. There- 
fore, these expressions become 



a 

1 + a 


e c la + 


1 

1 + a 




1 + a 




1 + a 


p< v 


(9) 


In order for this simplification to be correct, it is again neces- 
sary to impose the restriction — > — - . 

c 1 + 2a 


The quantities 

lated for values of 
table II for several 


pi v p< v 


f)c Zp , and 


have been calcu- 


0 = y = 1 in steps of 0.1 and are presented in 

if 


values of pz/c. 
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Loads due to aileron deflection . - The section coefficients of lift 
and twisting moment. 


and 


( pC l) a " 


( PC *)a 



( 10 ) 


respectively, are found for constant-chord trailing-edge ailerons of 
arbitrary length. All gaps between the wing and aileron are considered 
to be sealed. The calculation of the aileron derivatives 0c^ and 


is very similar to that of the indicial structural loads and is 


pc m 5 

included in appendix B. For most reasonable aileron configurations, if 
the aircraft is flying at supersonic speeds sufficiently high to satisfy 

1 > 


the Mach number limitation previously imposed ^that is, 


\ c 1 + 2a 

an aileron deflection on the left wing produces no loads on the right 
wing. For this reason, only the case is considered wherein the loads 
on one wing are independent of the aileron deflection on the other. In 
this case, the limitation on the modified aspect-ratio parameter pz/c 
for the analysis of aileron loadings in appendix B is that pz/c must 

c /c 

be greater than — — . Numerical values are given in table HI 


1 + 


2a -M* 

Z 


for 0 ^ ^ 1 in intervals of 0.1 for several values of pz/cj the 

computations have been made for b a /z = 1.0 and c a /c = 0.2. 


Significance of Mach number limitations . - Although the restrictions 
that have been placed upon the modified aspect-ratio parameter limit the 
utility of the aerodynamic results contained herein, these restrictions, 
in reality, are not serious. Ibis fact is substantiated by considering 
a typical configuration such as that used for the example contained in a 
subsequent section. For this wing, which has an aspect ratio of 3*6, 
the pertinent parameters are z/c = 1.5* a = 0.2, c a /c = 0.2, and 

b a /z = 1. The restriction imposed upon the expressions for the loads 


82 X 

due to structural deformation and roll, — > , and that for the 

' c 1 + 2a' 

loads caused by aileron deflection, — > — - — , become M > 1. 108 

c b= 

1 + 2a - - 2 . 

L 

and M > 1.05^-, respectively. Since the validity of linear aerodynamic 
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theory is questionable near a Mach number of 1, these limitations are 
of little consequence. 


Aeroelastic Solution 


Structural and rolling -moment equations .- In the solution of the 
aeroelastic rol lin g problem, not only must structural equilibrium 
(eq. (l)) be satisfied but also the equilibrium of moments about the 
rolling axis. If the loads in equation (l) are replaced by the loads 
arising from the various causes as derived in the preceding section, 
the equation specifying structural equilibrium becomes 


e(y) 





de 


d£ + 


c GM(y,Ti) jT pc^ObO H d£ 

y f^/ o p % (n) + c 


dTj + 


P c z 5 M + c GmCy^) PCuigCn) 


dq (11) 


Rolling-moment equilibrium is attained by setting the total rolling 
moment equal to zero. This condition can be written as 

a 


f (a Z + T]) L(tj) dll = 0 
J 0 


Again, the results of the preceding section can be vised to give 


qc 

0 


f 0 (aZ + ^ f Q p % (lb ^ + §7 P c z p (n) + 5 P%(n) 


dr) = 0 


( 12 ) 
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Alternate structural equation . - Equations (ll) and (12) completely 
express the necessary conditions for this aeroelastic problem. It is 
to be noted that the quantities 5 and qc/p are, in general, known, 
the structural twist 6(y) is unknown, and the aeroelastic rolling 
rate pb/2 V is the quantity that is desired. The simultaneous solu- 
tion of these equations in closed form is, in actuality, extremely 
difficult if not impossible to obtain. Some sort of numerical solution 
is therefore inevitable. One type of numerical solution, which utilizes 
a collocation technique in the solution of equations (ll) and (12), is 
derived in a subsequent section. 

The numerical methods with which the indicated integrations in 
equations (ll) and (12) are performed can be expected to be accurate 
when applied to reasonably well-behaved functions. In addition to 
numerical integration, however, numerical differentiation must be used 
in order to express the rate of twist d0/dy, upon which the loads due 
to structural deformation depend, in terms of the twist 0. As is well 
known, the process of numerical differentiation is not as accurate as 
that of numerical integration and therefore should be avoided if possible. 
The numerical differentiation can be eliminated by differentiating equa- 
tion (ll) with respect to y, which operation yields 


e=j („,e) § at + 

dy P J Q dy J Q iQ d £ 


dGM(y,T]) 

C dy 

qc pb r l 
y 2Vj n 


I 0 2 ^ 


dt| + 


dGj/y,!]) dG^y,^) 

+ 0 —5; — 


dy 


dr) + 


qc 


f 

J 0 


dGL(y,n) „ . dG M (y,n) 

P % (,1) + c — pc V n) 


dy 


dr] 


(13) 


In this manner, equilibrium equation (ll) is replaced by equation (13) 
and the twist 0 no longer appears in the problem. Hie rate of twist 
becomes the unknown and can be used directly in the numerical solution. 

The use of equations (12) and (13) is particularly desirable when 


the functions 


dG L (y,n) 

dy 


dG^y,!}) 

dy 


, which are merely the rates of 


and 
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twist due to concentrated loads and torques, respectively, can "be found 
accurately. The example considered in a subsequent section is one for 
which it was possible to obtain these functions analytically. When 
complicated structures are to be dealt with, however, recourse must 
usually be made to experiment or to approximate theories in obtaining 
the structural properties. For these structures, probably only discrete 
influence coefficients for twist can be found; the deter mina tion of the 
rate-of -twist influence functions would be subject to at least the same 
inherent inaccuracies as the numerical differentiation discussed pre- 
viously. In such cases, therefore, the use of equation (13) instead of 
equation (ll) would not be advantageous. 

Aileron reversal . - If the rolling-moment equation (12) is solved 
for 5, the result is substituted into equation (ll) or equation (13), 
and pb/2 V is set equal to zero, the following equations result: 

For equation (ll). 
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For equation (13) , 



Either of these homogeneous integral equations expresses the condition 
for aileron reversal, where the dynamic pressure at reversal q rev - 

appears as an eigenvalue. 

Rol 11 ng effectiveness . - The effect of the elasticity of the wing 
on the rolling behavior may he determined by examining the ratio between 
the rate of roll of the flexible wing and the rate of roll which would 
occur if the wing were rigid. This ratio, the "rolling effectiveness," 
is designated as f>, where 

(pb/2V) F 

? (pV 2V )r 

In this equation (pb/2V) F is the rate of roll determined directly 
from the aeroelastic equations. The quantity (pb/2V)^, the rate of roll 
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for the rigid wing, is found by setting the twist per unit length equal 
to zero in equation (12): 



f (az + y) pc, (y) dy 



J ( a * + y) 0c, (y) dy 


(16) 


Simplifications for wings with an elastic axis .- If the wing is 
constructed in such a manner that an elastic axis does exist as previ- 
ously described, the expression for structural equilibrium is given by 
equation (2) instead of equation (l). It is to be remembered that, in 
this case, the twists are dependent only on Q(y) , the torque about the 
elastic axis, and not on the load L(y) . The aeroelastic structural - 
equilibrium equation which results from using equation (2) instead of 
equation (l) can be obtained by merely deleting the terms in equa- 
tion (ll) that involve G-^(y , T]) and replacing the section moment coef- 
ficients about the midchord by the section moment coefficients about the 
elastic axis. Thus, the derivatives pc^ , and pc^ should be 


replaced by 



pc n , and 



respectively, where 


P c q 0 ~ P'thq “ e P c z Q 
0eq p = ficmp- e 0c Zp ^ 


(17) 


In these equations, as before, e is the distance measured forward from 
the midchord to the elastic axis, expressed as a fraction of the chord. 

The equations analogous to equations (13), (l4), and (15) are 
obtained in exactly the same manner. 


NUMERICAL METHOD 


Since an exact solution of the aeroelastic equations is not feasible, 
even for the simplest configuration, this section is devoted to the 
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presentation of a numerical method of solution. The method is based on 
satisfying the equations at a number of discrete spanwise stations. 
Matrix notation is used as an aid in organizing the numerical procedure. 
A step-by-step outline of the procedure is included at the end of this 
section in order to aid the reader. 

The accuracy of the numerical method is directly dependent upon the 
number of stations used. Experience has shown that for this particular 
problem 11 stations, defined by the end points of 10 equal spanwi.se 
intervals, are sufficient. The derivation that follows is therefore 
based on 10 equal intervals, the extension to other numbers of intervals 
being evident. Simpson's rule is used to perform the integrations and 
parabolic dif ference equivalents are used to replace any necessary 
derivatives . 


Matrix Operations 

Two distinct types of integration appear in the aeroelastic equa- 
tion. The first is of the form 



and the second is of the form 



dT] 


The integrands in these equations are evaluated at each of the spanwise 
stations t) 0 , t^, t\ 2 , ... r\ ± , .. . nio> where, for equal intervals of 

width e, m = ie. In a s imilar manner this subscript notation is used 

to denote evaluation of the integrands at each spanwise station. Thus, 
f i = f M = f ( ie )> g i = s(n ± ) = g(ie), and h^ = = h(ie,je). 

For 10 intervals, e = l/lO and Simpson's rule becomes, for each of the 
integrals, respectively. 



in = ^( f 0«0 + ^ f lSl + ^2^2 + 


^383 + . . • + ^fggg + f 10610) 
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and (for any single y-jJ 


f 0 H(y ± ,il) g(n) dn = ^(b 10 g 0 + 


+ . . • + ^iggg + h i,10 g lo) 


In matrix form these integrals can he written as, respectively. 


0 


f(n) g(n) in = £i]|Xll g i 


and (for all y^) 


f h(y,T]) g(n) in 

Wo i 


*±J 


g J 


The row matrix is given hy 


|*i “ f 0 f l f 2 f 3 • ' * f 9 f 10 


g l 

g 2 

Sx 


g i 




g 


10 


and a typical column matrix by 



18 


NACA TO 3067 


where the subscripts denote the position of each element. The square 
matrix is given by 








^00 

tod 

too2 

••• too, 10 


h 10 

ton 

to]_2 

top , 10 

II 

1 -ral 

u£j 

^20 

• 

*21 

• 

CVI 

^ * 

tog;, 10 

• • • • 


h 10,0 

toio,l 

h 10,2 

topo, 10 








where the first subscript designates the row and the second subscript 
designates the column in which an element appears. This subscript nota- 
tion is used throughout this section. The integrating matrix is given 
toy 



(18) 


All derivatives which appear in the aeroelastic equation are of 
first order; difference equivalents based on passing a parabola through 
three adjacent points are used herein to approximate the derivatives. 
For the points i = 1, 2 , ... 9 the standard difference equivalent 
derived by finding the slope at the center of the three ordinates is 
used: 



2e 
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For the end points i = 0 and i = 10, the slope at -the exterior of 
the three ordinates must he used: 


df\ = -5 f 0 + 1 ~ f 2 

W 0 2e 



f 8 " ^ f 9 + 5f 10 
2e 


In matrix form with e = l/lO, these expressions become 



where the differentiating matrix 



-3 4 -X 
-10 1 
-10 1 
-10 1 
-10 1 
-10 1 
-10 1 
-10 1 
-10 1 
-10 1 
1-4 3 


(19) 


Aeroelastic Matrix Equation 

The aeroelastic equations may now he written in matrix form by 
using the results of the foregoing numerical analysis of integrating 
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and differentiating processes. Equations (ll) and (12) "become, 
respectively. 



( 21 ) 
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If these equations are divided through by 5 and the resulting 
matrix equations are combined, the following partitioned matrix equa- 
tion results: 


- 1 ~ 

-kL 1 ! -Mil* 
T ! 


| 0 k/s| 


Ci 

LM i E 
_ 1 _ 


'pJoj 


F 


( 22 ) 


where 



is the unit matrix. 


The submatrices are given by 
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where 


E = - |aZ + yjJ 

S 

( ?Cj p) 3 

F = aZ + yj 

S 

O 30 ^ 




i, j, k = 0, 1, 2, . . . 10 


For a given configuration flying at a given speed and altitude, 
all the quantities in equation (22) are known except the twist 0 at 

each station and the rate of roll ~ -^8; solution of this matrix equa- 
tion yields these quantities. 

The rigid rate of roll, as found hy expressing the integrals in 
equation (l6) in numerical form, is merely 

m s - 1 


(23) 


If the rate -of -twist influence coefficients are known, it is prefer- 
able to employ equation ( 13 ) instead of equation (ll). In this case, the 
partitioned matrix equation analogous to equation (22) is 


1 

qc 

P 

1 1 
H 

1 l 

>i; 

®i 


/de\ 

W/k 

5 

= 

— 

O 

I- 1 - 

— 


% 


E 


2V/ 


F 


(210 


where the submatrices are given by 
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and the expressions for the scalars E and F remain the same as those 
given immediately after equation (22). 

By applying the same numerical processes, the aileron-reversal equa- 
tion (l4) (or eq. (15)) can he put in matrix form. Equation (14), in 
which 0 is the variable, becomes 



( 25 ) 


The submatrices are defined immediately after equation (22) . 
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Similarly, equation ( 15 ) , where d0/dy 


is the variable, becomes 



( 26 ) 


The submatrices are defined Immediately after equation (24). 

A solution to equation ( 25 ) (or eq. (26)) is easily found by using 
matrix iteration. The process converges to give the lowest eigenvalue 
from which the dynamic pressure at aileron reversal may be obtained. 


Eeduced Matrix Equations 

Up to this point, the numerical analysis has been based on the use 
of 11 spanwise stations, a number that was deemed necessary in order to 
obtain the desired accuracy because of the ill-behaved nature of some 
of the aerodynamic loading functions . In most cases, however, the twist 
(or rate of twist) is well behaved and, therefore, requires fewer sta- 
tions for adequate specification. Utilization of this fact allows a 
considerable saving in the amount of work necessary to solve the matrix 
equations because of the fewer degrees of freedom involved. 

If either the twist or the rate of twist is specified at the even- 
numbered stations, an interpolation procedure can be used to determine 
the values at all 11 stations. The particular type of interpolation 
used herein is obtained by passing a fifth-degree polynomial through 
the even-numbered stations and then evaluating this polynomial at the 

This interpolating procedure can be written in 

0k I = M 1 91 

si 

. 10 and Z = 0, 2, 4, •. . . 10 and the inter - 
is given by 



odd-numbered stations, 
matrix form as 


or 


where k = 0, 1, 2, . . 
pointing matrix N 
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1 

0 

0 

0 

0 

0 

0.246094 

1.230468 

-0.820312 

0.492187 

-0.175781 

0. 027344 

0 

1 

0 

0 

0 

0 

-0.027344 

0.410156 

0. 820312 

-0. 273437 

0. 082031 

-0.011718 

0 

0 

1 

0 

0 

0 

0.011718 

-0.097656 

0.585958 

0.585958 

-0.097656 

0. 011718 

0 

0 

0 

1 

0 

0 

-0.011718 

0.082031 

-0.273437 

0. 820312 

0.410156 

-0 . 027344 

0 

0 

0 

0 

1 

0 

0.027344 

-0.175781 

0.492187 

-0. 820312 

1.230468 

0. 246094 

0 

0 

0 

0 

0 

1 


The twist (or rate of twist) now needs to he known only at the 
even -numbered stations , the quantities at the odd-numbered stations 
being obtained by interpolation. Therefore, only the even rows of the 
influence -function matrices are necessary; these matrices become rec- 
tangular with 11 columns and 6 rows. When these simplifications are 
used to write the matrix equations for equations (ll) and (IP) and the 
resulting equations are combined, the following partitioned matrix equa- 
tion results: 






V 




(28) 
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where the submatrices are given by 



Note that the subscripts j and k assume all values, both even and 
odd, and the subscripts m and l assume only the even values. 
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The analogous partitioned matrix equation to he used when the rate- 
of -twist influence functions are known is 



where the suhmatrices are given hy 
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The sc alar s E and F are given after equation (28). Here again the 
subscripts j and k take on all values, whereas m and l assume 
only even values. 

The aileron-reversal equations may also "be expressed in this 
reduced-order form. Thus, equation (14), in matrix form, is 



^rev c 

P 





9 l 


(30) 


The submatrices are defined after equation (28). Similarly, equation (15) 
"becomes 



The submatrices in this case are given immediately after equation (29) • 
Again, standard matrix iteration procedure may be used to solve equa- 
tion (30) or (31) , tie amount of work "being approximately one-fourth as 
much as is involved in the iteration of equation (25) or (26) . 


Simplifications for Wings With an Elastic Axis 


Two operations are required to modify the aeroelastic matrix equa- 
tions when an elastic axis exists. First, all terms involving the 


matrices 


°L 


or 


dy 


are deleted. Second, the section-moment- 


coefficient matrices 


Pc 


*0 


PCr 


, and 


Pc 


% 




, and 


Sc ”s 


are replaced "by 


in accordance with equation (17). 
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Computational Procedure 

As lias teen previously mentioned, the task of calculating the 
rolling effectiveness for the given configuration flying at a given 
altitude and speed is straightforward; only the solution of a set of 
simultaneous equations is necessary. If a complete knowledge of the 
aeroelastic behavior over large ranges of speed and altitude is desired, 
however, the manner in which the various parameters are involved should 
he taken into account in determining the sequence of calculations. It 
should he noted that the Mach number enters the problem in a complicated 
manner, whereas the dependence on altitude is rather single; that is, 
the Mach number affects the aerodynamic matrices and the altitude affects 
only the dynamic pressure q. It is therefore obvious that the most 
economical way to perform the computations is to calculate the variation 
of rolling effectiveness with dynamic pressure for several constant values 
of Mach number. In addition, since the roll in g effectiveness is of little 
interest when the controls are reversed, the range of dynamic pressures 
should be restricted to values less than q rev . He dynamic pressure at 

reversal should, therefore, be determined for each value of Mach number 
before proceeding with the calculation of rolling effectiveness. 


An outline of the steps required in the determination of the aero- 
elastic effect on roll for a range of altitude and Mach number is 
included herein. For simplicity, only one of the several numerical 
approaches derived in this section - that is, the one wherein the twist 
influence functions G^ and and the interpolation procedure are 

utilized - is illustrated; the others follow the same outline, differing 
only in detail. 

(l) Evaluate, either analytically or experimentally, the twist 
influence functions Gp(y,Ti) and Gj l j(y, , n) at stations y/z = 0, 0.2, 

0.4-, . . . 1.0 due to loads and torques applied at stations q/z = 0, 
0.1, 0.2, . . . 1.0. From these values, form the ll-by-6 matrices 

[Nj “ 4 [Mnj]- 


(2) For a given Mach number (one should be chosen which results in 
a value of 0z/c which appears in the tables) determine the indicial 
section lift and moment coefficients 3°Zg Hese coeffi- 

cients, which are found by applying equation ( 7 ) to the values of P c Zq' 


and 


in table I, can then be used to form the square matrices 


( pc Z 


0 /jk 


and 


( pCm e) 
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(3) In a s imilar manner, form the column matrices 


( P0 *p) 


pc 


( PCm p); 

7 and 8c-, 
l V Q 2 a, 


and 


, making use of equation (9) in conjunction with the values of 


in table II. 


(4) Determine the section loading coefficients due to a unit aileron 
deflection. Use these coefficients to construct the column matrices 


( p %) 


j 


and 


( PCffl 5)- 


(These coefficients are listed for one 


aileron configuration in table III. Equations for the coefficients for 
a rather general aileron configuration are included in appendix B. ) 

(5) Compute the matrices [ w] , 1^*1' I 0 **]' and [_vj and 

the scalars E and F defined inmediately after equation (28) . The 
integrating and differentiating matrices |^S^j and Qd J in these defi- 
nitions are given by equations (l8) and (19); the interpolating matrix 
i- s given by equation (27); the row matrix jjiZ + y^J 1 b made up 
of the moment arms about the rolling axis. 


(6) Form the matrix 


Am* 


+ i 

F 


Cm* 


V 


(7) Obtain q rev by iterating this matrix (see eq. (30)). 

(8) For each of several values of q between zero and , form 

the matrix equation (28). Solution of this equation yields 

for each value of q. Obtain the rolling effectiveness 0 by dividing 

(§/i * (i/i ((§A) r - §)• 

(9) Repeat steps (2) to (8) for several other values of Mach number. 


SAMPLE APPLICATION 


In this section, aileron effectiveness and reversal speed are found 
for a specific aircraft by the approach set forth in this paper, which is 
termed "lifting surface theory," and by two simplified methods. Each half 
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of the exposed wing consists of a uniformly thick aluminum plate with 
an aspect ratio z/c of 3/2 and a thickness ratio t/c of 0.02. The 
full-span aileron is formed by bending the plate along the 0.8-chord 
line} therefore b a /z = 1.0 and c a /c = 0.2. The wings are mounted 
rigidly on the body, which has a radius of one-fifth the exposed wing 
semispan} therefore a = 0.2. These ratios are sufficient to define 
completely the configuration, since, as is seen later, the results are 
independent of absolute dimensions. 

Although the plate is considered to be bent in order to form the 
aileron, the plate is assumed to behave structurally as if no bend had 
been made. On the basis of this assumption the necessary structural 
influence coefficients are determined in appendix A by means of a sim- 
plified flat-plate theory. In this case, as could be expected, not only 
does an elastic axis exist but also this elastic axis coincides with the 
midchord. 


Calculation by Lifting-Surface Theory 


Since the influence functions are obtained analytically, the rate- 
of -twist influence function rather than the twist influence function has 
been determined in order that no numerical differentiation be necessary. 
In addition, the wing structure being free of discontinuities, the inter- 
polation procedure developed in the preceding section may be successfully 
employed. Equations (2 9) and (31) , modified as outlined in the preceding 

section, are therefore .used to compute and l rev , respectively. 


Before proceeding with the solution of equations (29) and (31) , the 
dimensional character of the structural, aerodynamic, integrating, and 
interpolating matrices involved in these equations should be considered. 
The aerodynamic and interpolating matrices are clearly dimensionless. 

-moment-arm matrix J^aZ + yjj 


The integrating matrix 


S 


and the rolling 


are proportional to the exposed semispan Z. The influence -function 
matrix for the example configuration, which is presented in table IV, 


written as the product of the quantity 



and a nondimensional 


is 


matrix. 


The aforementioned independence of absolute dimensions can now be 
demonstrated. Equations (29) and (31) can be written in the form 
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where the multiplication of the various submatrices by the indicated 
quantities yields nondimens ional results. Inspection of equations (32) 
and (33) shows that only ratios of the dimensions are involved, and that 
the quantity 


qcl 2 



is an important nondimens ional parameter for this particular configuration. 

The step-by-step procedure previously outlined has been followed 
for the example configuration and the results are included in figure 3* 
Since the dynamic pressure in this case is essentially an altitude param- 
eter, the results in figure 3 are shown plotted against the pressure ratio 
P h /P Q , where Pq is the standard pressure at sea level, as obtained from 

reference and P^ is the free-stream static pressure at altitude and 
is related to q by the equation 


■?h = 


2q 

7M 2 


in which 7 is the ratio of specific heats of air. 
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Simplified Methods 

Two simplified methods of analysis have also teen used to solve 
for q rev and 0 for the example. They differ from the method pre- 
viously described only in that the airloads are determined by simplified 
means. In the first method, lifting-surface theory is modified in a 
manner similar to that used in reference 6, wherein it is assumed that 
the chordwise center of pressure due to structural deformation and roll 
coincides with the elastic axis of the wing. Thus, the total twisting 
mQiuent about the elastic axis is equal to the twisting moment due to 
only the aileron, and equation ( 15 ) becomes 


de = qcf 6 r z 5G M(y^) 

<3y 0 Jq ^y 


0^(1]) dn 


(3*0 


The rolling-moment equation, equation (12), is unchanged. A solution 
to the problem has been obtained by using matrices in a manner similar 
to that employed previously. In this case, however, de/dy is given 

explicitly by equation (3*0 and, therefore, ^^5 can be calculated 

directly, the solution of a set of simultaneous equations being 
unnecessary. 

In the second method of analysis, the aerodynamic terms are derived 
on the basis of two-dimensional (strip) theory. For the case wherein 
the elastic axis lies on the midchord, the center of pressure due to 
structural deformation and roll coincides with the elastic axis of the 
wing and the only twisting moment about the elastic axis is that produced 
by the aileron. Therefore, equation (3*0 expresses de/dy exactly, and 
the rolling rate obtained by solving equations (12) and (3*0 is exact. 
Because of the simplicity of the strip theory, these calculations can be 
performed analytically. 

The results obtained by using these methods, termed, respectively, 
"modified lifting-surface theory" and "strip theory," are also shown in 
figure 3 for comparison with the results obtained by using lifting-surface 
theory. 


RESULTS AMD DISCUSSION 


The results shown in figure 3 indicate that the variation of rolling 
effectiveness 0 with the pressure ratio P^/Pq is> f° r practical pur- 
poses, linear. This linearity suggests the possibility that, for configu- 
ations of the type considered, only the computation of pressure ratio at 
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reversal (P^/Po) would "be required; the rolling effectiveness for 

smaller values of P^/Pq could, then he obtained by linear interpolation. 

Tn this way the computation time could be greatly reduced, the calcu- 
lation of (Ph/ p o) r ev ^ eing a relatively simple process. Some care 

should be taken , however, in making use of this linearity, since it 
appears to depend on the proximity of the elastic axis to the chordwise 
center of pressure due to angle of attack. For wings wherein the elastic 
axis is distant from this center of pressure, considerable curvature of 
the rolling-effectiveness curves could result. At the other extreme, 
the results in figure 3 obtained by using the two simplified methods are 
exactly lin ear, since, in both cases, the elastic axis and the center of 
pressure are coincident. 

An additional consequence of this linearity is that the accuracy 
with which a particular method predicts (Ph/-^o) re v a < ^- rec 4 measure 

of its ability to predict rolling effectiveness. A comparison of the 
values of (^hl^o) rev obtained by the two simplified methods with those 

obtained by the method presented in this paper is therefore shown in fig- 
ure 4, wherein (Ph/ p o) re v is Mach number. Although the 

results show very little difference in the values of (Ph/^o) re v aS 

obtained by the three methods at high values of M, considerable differ- 
ence exists at low values of M. The results obtained by the use of the 
modified lifting-surface theory are consistently unconservative; the 
results obtained by the use of strip theory are consistently conserva- 
tive. Modified lifting-surface theory neglects the twisting moments 
arising from the twist of the wing; only the torques caused by aileron 
deflection are considered. Neglect of the twisting moments due to angle - 
of -attack changes evidently reduces the resultant angle of twist and 
therefore reduces the adverse ro llin g moments caused by the deformation. 
When strip theory is used, the elastic twist is again lower than that 
obtained by lifting -surface theory because, again, only the twisting 
moments caused by the aileron deflection are present. Here, however, 
the absence of the finite -span effects actually results in a greater 
adverse rolling moment, even though the twist causing the adverse rolling 
moment is in itself smal l er. 


Further mention should be made of the behavior at large Mach numbers. 
Actually all three methods should yield the same results as M approaches 
infinity because, as M increases indefinitely, lifting-surface theory 
approaches strip theory. As an illustration of this fact, the rate of 


roll for a rigid wing 



is shown as a function of M in figure 5* 



R 


The values of 


obtained by using lifting-surface theory rapidly 
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approach the values obtained hy using strip theory as the Mach number 
increases. Thus, at high values of M the use of strip theory should 
yield aeroelastic solutions of an accuracy comparable to that obtained 
by using lifting-surface theory. 

Tn actual application it may be convenient to have rolling effec- 
tiveness given as a plot of pressure ratio against Mach number for con- 
stant values of rolling effectiveness. A graph of this type may be made 
by cross -plotting the information contained in figure 3 and is included 
in figure 6. Also included in figure 6 is an additional ordinate that 
gives the altitude as obtained from the standard-atmosphere table in 
reference 5* 

The rolling effectiveness 0 at any time during a particula r flight 
may be determined as a function of Mach number if a history of the flight 
is known in the form of a plot of standard altitude against Mach number. 
For example, consider two constant-altitude flights, one at 30,000 feet 
and the other at 20,000 feet. Hie resulting variation of the rolling 
effectiveness of the example configuration with Mach number is shown in 
figure 7* 

The indicial-solution (aerodynamic-influence-function) approach used 
herein in the calculation of the aerodynamic loads exhibits considerable 
promise of being applicable to plan forms other than rectangular. The 
delta wing and the law-aspect-ratio swept wing at supersonic speeds can 
be handled in the same way as the rectangular wihgj care should be taken 
in these cases, however, to account for chordwise bending where necessary. 
Other static aeroelastic problems such as center-of -pressure shift (which 
has been considered by Frick and Chubb in ref. 1 far high-aspect-ratio 
swept wings) and load distribution seem to be amenable to analysis by the 
methods contained herein. It might even be possible to extend the 
approach to take into account unsteady aerodynamic effects and thereby 
to obtain accurate solutions to flutter problems. The calculations in 
the last case would undoubtedly be arduous and the main usefulness of 
the approach would be to establish a basis for the evaluation of mare 
practical but necessarily less accurate solutions of the flutter problem. 

Problems involving configurations about which the flow is not sub- 
stantially potential are generally not amenable to this type of approach. 
More specific ally , the success of the approach depends on the applica- 
bility of linearization to the aerodynamics and on the ability to calcu- 
late the necessary indicial load distributions. 

COLLUDING REMARKS 

A method has been developed for the prediction of the aeroelastic 
effects on the roll of rectangular wings in supersonic flow. The method 
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is 'based on the use of influence functions, either analytically or 
experimentally obtained, to calculate the structural distortion due to 
airloads. The airloads, themselves, are calculated on the basis of line- 
arized lifting-surface theory by superposing basic airloads resulting 
from elementary angle -of -attack distributions. The solution of the aero- 
elastic equations has been obtained by means of a numerical procedure 
suitable for use with desk- type calculators. 

Results for an example configuration indicate that the variation of 
rolling effectiveness with the free-stream static pressure at a constant 
Mach number is almost linear; a good approximation may be made by assuming 
a linear variation. Thus, in any other cases wherein this linearity can 
be expected - that is, when the elastic axis is near the center of pres- 
sure due to angle of attack - the calculations may be greatly simplified 
in that it is necessary to compute only the free-stream static pressure 
at aileron reversal. 

The results obtained by using the method of analysis presented in 
this paper are compared with the results obtained by using simplified 
methods of analysis. Although aerodynamic strip theory is valid at high 
Mach numbers, too conservative results are obtained at low supersonic 
Mach numbers because of neglect of finite-span effects. A modified 
lifting-surface theory, in which twisting moments due to structural 
deformation and roll have been neglected, yields results which are 
unconservative at low supersonic Mach numbers. 


Langley Aeronautical Laboratory, 

Rational Advisory Committee for Aeronautics, 
Langley Field, Va. , December 2, 1955 • 
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APPENDIX A 

STRUCTURAL ANALYSIS OF THE EXAMPLE CONFIGURATION 


The deformations of the wing of the example configuration are ana- 
lyzed herein "by using the approximate plate theory presented in refer- 
ences ^ and 7* This method expresses the deflection w of a plate at 
any point (x,y) in the form w(x,y) = W(y) - x 0(y) where W(y) is the 
deflection of the wing at the midchord and 9(y) is the twist of the 
wing. (See fig. 1.) An energy solution of the problem is used wherein 
the expression for the potential energy of the plate is written in terms 
of the approximate deflection function. Minimization of this potential- 
energy expression yields two ordinary differential equations in W(y) 
and e(y) . Since the deflection W(y) has no effect on the airload, it 
is eliminated from the two equations; the process yields a single equa- 
tion in 0(y). This equation has been derived in reference 7 (eq. (A22) 
of that paper) and for a rectangular cantilever plate of constant thick- 
ness t becomes, in the notation of the present paper. 


— 2^2 _ 2(1 - p)Dc 22 
12 dy 3 dy 



M(5) d£ 


(Al) 


with the boundary conditions 


0 ( 0 ) = 


d0(o) d 0(l) 


<3y 


= 0 


dy £ 


In equation (Al), D is the plate stiffness: 


D 


EV 


12(l - p 2 ) 


where E is Young's modulus and p is Poisson's ratio. As could be 
expected for this structure, the twists are seen to be dependent solely 
on the twisting moment about the midchord M(y) . 

The solution to equation (Al) for a general M(y) can be obtained 
by superposing indicial solutions found by considering the moment to be 
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concentrated at any spanwlse station q. Thus, let 

M(y) = S(y - n) 

where 5(y - q) is the Dirac delta function. 

Equation (Al) now he comes 

b£ ^. 1 ) _ 2(1 , )a> ^ . . i( „ . y) - m 

12 ^ 


where _^(q - y) is the unit step function: 



I 

II 

H 

(y £ n) 


1 

II 

O 

(y > n) 

fphp rate-of -twist influence function =r is the quantity required 

oy 

in order to obtain the rate of twist by superposition. This quantity is 

SGuCy^ti) 

1 L. h > _ ™„ nfl - r 

-•> 

) - 


8y 

Dc 5 A 2 cosh AL 




sinh A^l - Y^sinh A 

(y ^ n) 

> (A3) 

8%(y,q) 

_ 121 ^ - sh 4 - % , , , 2 ) 

(y > 11 ) 


8y 

Dc 5 A 2 cosh A 

, 


where 


A = i\|24(l - p) 
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With regard to this solution, it should be noted that the quantity 
12 Z^/dc^A^ can be written in the form 1 w ^ere G is the shear 
modulus of elasticity of the material. 

Superposition of the rate-of -twist influence function yields 


Mk) . r 1 M(n) d» 

dy Jq By 


(Ah-) 


Values of BGy/dy for z/e - 3/2 and p = l/3 (which yields 
A = 6) have been computed for 0 ^ y ^ 1 in intervals of 0.2 and 

0 ^ ^ 1 in intervals of 0. 1} the results are included in matrix form 

in table IV. 
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APPENDIX B 

DERIVATION OF AERODYNAMIC COEFFICIENTS 


In this appendix the aerodynamic loads necessary for computing 
rolling effectiveness are derived. The loads are found in the form of 
section lift and moment coefficients hy applying linearized supersonic 
lifting-surface theory to the rectangular wing with three different 
ahgle-of -attack distributions: the unit-step angle -of -attack distri- 

bution shown in figure 2(c) , the angle -of -attack distribution caused by 
ro lling , and the angle -of -attack distribution resulting from aileron 
deflection. 

Analyses of each of these problems are contained in the literature. 
The unit-step problem is essentially the same as the problem of finding 
the loads due to deflection of an outboard aileron, which has already 
been solved (see, for example, ref. 8)$ the rolling problem has been 
treated by many investigators (see, for example, ref. 9 )* There does 
not seem to be any report, however, that gives the desired coefficients 
in a form sufficiently complete for the purposes of this paper. For 
this reason and also for convenience, the necessary aerodynamic quan- 
tities are derived completely herein. 

When linearized lifting-surface theory is used, the lift per unit 
area of a thin wing is given by 


p(x,y) 


V 8x 


(Bl) 


where x is positive in the direction of the airstream. The poten- 
tial cp, evaluated at the surface of the wing, is 


<P(*,y) 


V 

it 



q(g,£) dg d£ 
l/(i - x ) 2 - p 2 (£ - y ) 2 


(B2) 


where cr(|,£) is the local angle of attack of the wing. In general, for 
a rectangular wing, the region of integration S includes the entire area 
on the wing within the forward Mach cone from the point (x,y) . However, 
in order to obtain the potential at a point near a wing tip, the proper 
region of integration S is determined by using Eward's method (ref.10). 
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in 


The spanwi.se loading and twisting moment ahout the midchord are 
given by 


L(y) 

M(y) 



(B3) 


If the expression for p(x,y) is substituted into equations (B3), and 
the results are integrated, the lift and moment become: 



where the fact that the potential is zero at the leading edge has been 
taken into account. 


The foregoing equations can be nondimensionalized by letting 


P C Z 


Lp 

qc 
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The resulting equations are 


p c z(yi) = it ®( 1 >yi) 


Pc m( y l) = “T 


| ®(^yi ) - / ®( x i*yi) a*: 

a Uq 


(B5) 


where 


®( x i>yi) = 


g (Su£i) d £i d li 

S \j(x 1 - ! x ) 2 - m 2 ( yi - 


(b6) 


and m is the modified aspect-ratio parameter pz/c. 


Derivation of P c Zg' an^- 3 c it1q ' 


The equations for the spanwise loading P c Zg' and the twisting 
moment about the midchord pc mg ' due to the unit-step angle -of -attack 
distribution shown in figure 2(c) become, from equations (B5), 


4m 


|k V( 3r i’ ,| i) * v 4, ( 1 ’ y i’ii) 


4-m 


^'(yi^i) =~ 


> 


(B7) 


\ ®(i,yi,Tii) - f ®(xi,yi,Tii) dx-L 
^ u 0 


where the potential $ is found for an angle -of -at tack distribution 
defined by the unit step function: 

tf(yi^i) = l(r\ 1 - yi) 


CT (yi^i) = 1 
a (yi^i) = 0 


or 


(41 ^ yi) 
K < yi) 
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Equation (b 6) becomes 

4 (xi,yi,Tii) = JJ - = £l = (b8) 

e ^(*i - ii ) 2 - m 2 (y x - £i ) 2 

The region of integration Z includes only that portion of S which 
lies within the region of the wing between the wing tip and the position 
of the unit step in angle of attack at tj-^. 


The plan form and angle -of -attack distribution considered are those 
shown in figure 2(c). From this figure it can he seen that the position 
of the unit step is restricted to the range 0 ^ ^ 1. If, in addi- 

tion, the modified aspect-ratio parameter m is restricted to values 
greater than — ^ — , the potentials 4> for the finite wing in fig- 

ure 2(c) are exactly the same as the potentials on the finite wing por- 
tion of the following semi -infinite plan form which has been obtained by 
allowing the left tip and the trailing edge of the finite wing to 
approach infinity: 



I 




' — ' 

s' 

TL 

\ 1 ^" 

DKaf H 


EC (a) ^ 

s' 


21(a) 

^ lb) 

— ^ 



x i 


Thus, attention is confined to the derivation of the potentials on this 
semi-ifafinite plan form. With the potentials known, the loading coeffi- 
cients on the finite wing are given by equations (B7) . 

The Mach lines shown in the above sketch bound a number of distinct 
regions. These regions are significant because the area of integration 
in equation (b 8) takes a different geometric form in eachj it is to be 
expected, therefore, that the potential $ is given by a different 
equation in each region. Hie proper area of integration Z (shown 
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shaded) and the potential 4 for the various regions, as found from 
equation (B8), are: 

In region I: 



4 


I 


£ 

m 


X 1 


In region II: 



®II = | ~ ^1 + x ! -tan" 1 ]j 


Wi 
x£ - 
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Notice that the only difference in the equations for $ju(a) an d 
$m(3) is that the radical in the inverse-hyperbolic-tangent term in 
the equation for $hi(a) reciprocal of the radical in the 

corresponding term in the equation for Thus, the following 

form holds for both regions: 


K m 


X3 tan 


-1 


I 


x i + m Pi - yi) 1 

— — + m/T)-, - y-^tahh - L k n 

x i - m ( T ii - yi) ^ 1 1} 1 


where 


k l _ 


x i - m (X. ~ W 
I x i + m pi - y±) 


(*i £ \) 


k i - 


x i + m 0 ii - yi) 


Pi - m 0i - yi) 

A similar situation occurs in the derivation of ^jv^a) 


(*i - ’ll) 


and $jv(B) 

and of ^V(a) and- ®V(B) * Therefore, the equations for and $ y 

may he derived for either the subregion (a) or (b) and applied through- 
out the entire region. 

In region IV: 



'IV 


iX, 


+ x i 


tan” 


i | x i + m (Bi - yi) 
x i - “(’ll - yi) 


tan - 


-1 Pi " "Vi 


myi 


+ m(r [ 1 - y-^) tanh -1 ]!^ 
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Since the spanwise loading pcj^' is determined by the value of 0 

at the trailing edge (see eq. (B7)), the equation for P c Zg' is 

dependent only on the region that the trailing edge is in. For example, 
if for a given value of y-]_ the trailing edge is in region III, the 

value of 0 c 2 g* is 

if ®m( y l' 1 ' ,1 l) 

Thus, the equations for $ c Iq' are: 



and the partic ular form to use for pc^ 1 is determined by the region 
in which the trailing edge falls. 
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^9 


The determination of pc^' requires evaluation of the integral 

of the potential along the chord of the wing. However, the indefinite 
integrals of the expressions for the potentials valid within each region 
happen to he continuous across the boundaries of each region. Further- 
more, the indefinite integral of the potential that is valid in the 
neighborhood of the leading edge is zero at the leading edge. Thus, the 
integral of the potential along the chord may be found by ev alua ti ng at 
the trailing edge the indefinite integral of the potential that is valid 
within the region in which the trailing edge is located. When the 
trailing edge falls successively in each of the regions I to VI, the 
resulting equations for pc' are, therefore. 
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where 


1 - m(n 1 - yi ) 
1 + m ( T l 1 - yi) 



(yi ^ ^1) 


1 + m(T] 1 - y^ 

- m (1i - y x ) 



( y i = ’ll) 


As "before, the particular form to use for 
region in which the trailing edge falls. 



is determined hy the 


Equations (B9) and (BIO) may now he used to obtain the spanwise 
loading and twisting moment about the midchord due to a unit-step angle- 
of -attack distribution for particular values of T]^ and m. In summa- 
rizing the information, it is convenient to separate m and tj 1 into 
various ranges. Sketches of the Mach lines and summaries of the par- 
ticular forms to use in each range follow: 


For m iT 2 


and 0 -g ri 1 ^ J: 
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Use the coefficients 


for region II when 0 % y, 1* — 

± m 

for region I when ^ = y^ = - i 

for region m when ~ ^ ^ ^1 ^ ^1 + 
for region VI when ’ll + m “ -^1 - ^(l + a) 


— - ^ m < 2 and 0 tu ^ i; 

1 + 2a “ 1 “ m 


/ \ / 

/ /\ 

/ m / rz ) 

/ / / 

/ / / 

/ / / v 

/ / / -X- 


-+T 1 - 

1 '1 m 


m ^1 


Use the coefficients 


for region V when ® - -^l = m ” ^1 

for region IV when — - ri. ^ y, ^ — 

0 m 'l^lm 


for region HI when — ^ y_ S — + ti_ 

m J 1 m ‘1 

for region VI when i + = y^ - 2(l + a) 


Bln 
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Use the coefficients 


for region H when 0 ^ y, ^ - -i 

1 '1 m 

for region IV when Ti n - i < y, £ i 

for region III when — ^ y_ < ri n + i 

m J 1 ‘1 m 

for region VI when + i % y 1 % 2(1 + a) 
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In tills section, the potential $ is derived for a finite rectan- 
gular wing with a nondimens ional span of 2(l + a) . In order to simplify 
the derivation, it is assumed that no point on the wing is influenced hy 

both wing tipsj this assumption is valid provided that m ^ — — — . 

1 + a 

The Mach lines originating at the intersections of the leading edge 
and the tips divide the wing into two distinct regions (see the following 
sketches) and the potential $ for these two regions is given hy equa- 
tion (B12) as follows: 


In region I: 



In region H: 
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*n- -- 1 


m 1 + a 


1 - 


yi 


- f + X i ( 1 - yi) te “- 1 |^^ yi 


I + -1 


The equations for the spanvise loading 0c z and the twisting moment 
about the midchord pc mp are found in a manner similar to tha t used to 
determine the equations for 0 c Z q » and Thus, 



The equations for obtained in a manner similar to that used pre 

vlously, are 


PS)r = ° 


0 %) 


_ 8 1 


n 


rt 1 + a 


K 1 - T + - -1 f 2 


- > (B14) 


8 a 
it 1 + a 


fj^i(i - “yi) 5/2 
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In order to conform to the assumption made in the aerodynamic 
analysis, the airloads over the rigid plate in the center section of 
the wing are to he neglected. Therefore, the limit of the modified 

aspect-ratio parameter can he lowered from m = — - — to m ^ — i . 

1 + a 1 + 2a 

Summaries of the particular forms to use and illustrative sketches 
follow: 

For m ^ 1: 



For i — = m ^ 1: 

1 + 2a 

y. 


Use the coefficients for region II when 0 ^ y^ ^ 1 
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It should, he noted that each expression for pc^ and pc^ can 
he separated into two parts - one with the coefficient ^ 


1 + a 


and the 


other with the coefficient 


- each of which is independent of a. 


1 + a 

The first part is the result that would he obtained on the right wing if 
the wing were rolling about the x-axis (see fig. 2(a)); it is therefore 
designated P c jp or pc^ • ® ie second part can he shown to he the 

section lift or moment coefficient which results from a unif orm angle 
of attack; it is therefore designated P c x a or Pc^. Thus, the quan- 
tities 0 Ct and PCmp can he written as 


P c Zp 


^ C 2-n = 

ip 


Pamp 


■T— = PC, 

1 + a ^p. 


1 + a 


P c 7 
t a 


TTT pCm P 0 " T+~& PCj % 


where the coefficients in the various regions are: 


- -M 1 - n> 
= 4 


(No), 


= 0 


Cx)j. ■ 0 


(N 0 ) 


n 


8 

x. 


1 - 


yi 


- ^ + i 1 - yi 




(p^a)^ - f + ^Jrr 

(^0)33: = ‘ M 1 ' T + - ^l) 372 

(p c %) n - - "Wi) 5/2 
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Derivation of pc^ and. P c mg 


Tn this section, the aerodynamic loads due to aileron deflection 
are derived for a trailing-edge aileron of constant chord c a and 
span b a . The spanwlse location of the inboard end of the aileron is 
left general; the outboard end of the aileron coincides with the wing 
tip. It is assumed that there are no gaps between the aileron and the 
wing. Furthermore, the portion of the wing forward of the hinge line 
may be ass ume d to be absent because the wing is considered to be at zero 
angle of attack, the aileron hinge line is unswept, and only supersonic 
speeds are considered. The resulting configuration is exactly the same 
as tha t encountered in the derivation of P c jg flCmg. In addition, 

since, for most reasonable aileron configurations, the pressures on the 
right wing are unaffected by the deflection of the left aileron, the 
total spanwise loading pc^g and the total twisting moment about the 

midchord pc^. can be found by properly substituting the aileron dimen- 
sions for the wing dimensions in the results already obtained for P c Zg' 
and pc m g 1 . In order that the foregoing be true , the restriction 
- c a-i 

m m is necessary. 

1 + 2a - b &1 


If c a is substituted for c and b &i = b a /l for r\ ± in the 


expressions for the loads due to structural deformation, the load due 
to a unit aileron deflection will result. Several more steps are neces- 
sary, however, in order to produce the desired results P c Zg an< i Pc^. 

In the first place, the reduced aspect-ratio parameter which was p l/c 
becomes pz/c a = m^c a ^, where c ap = c a/ c * ^- ri sedition to the substi- 
tutions mentioned in the preceding paragraph, therefore, the substitu- 
tion of m/c a -, for m must be made, both in the expressions for P c jg' 


aDd P C JHg ' 


find in the expressions which define the l im its of applica- 


bility of these terms. Furthermore, in the nondimensionalization of the 
moment expression the q uant ity c is involved. In order to preserve the 
nondimensionalizing coefficient qc/p for the aileron loads and qc^/p 
for the aileron moments, the expressions for P c lg ' (yi^a^) an< ^ 

“ust te multiplied by c^ and c ai 2 , respectively. 


Lastly, since the moments are taken about the midchord in the case of 
pressures due to the structural deformation, direct substitution yields 
the moment about the aileron midchord; thus, the moment must be trans- 
ferred back to the wing midchord. 
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The proper forms to use in determining pCj and (3^ are 


summarized as follows for values of m greater than "both 


' a l 


2c ; 


a l 


1 + 2a - h 


a l 


^an 

For h a < 1 ±: 

a l - m 


Ja l 


and 




1 

3ZI 

/ 

/ 

/ 

/ 

/ 

\ x / 
\ I / 
\ / 

m\ /n 


— Yzr 


'<*1 m 


m 


k _ 

Dq i m 


< < :*i 


for region II when 0 = y ^ 


m 


for region I when 


l l< , r < 


m 


s yi a ^ - 


c a. 


m 


Use the coefficients < 


a l 

for region HI when t a „ ^ y n < + 

1 m 1 1 

C a 1 

for region VT when h a ^ H — — ^ ^ 1 
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^ ^ C a 1 

for region II when 0 = y, = 

± m 


I Cg^-. 

Use the coefficients < for region I when ^ y x = h &i — 


for region HI when h a = y-. = 1 

1 m 


The expressions for the coefficients in the various ranges are: 
= 4c ai 
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TABIE I.- PARTIAL SECTION LIFT AND MOMENT COEF FICIENT S SCj e ' AND 0Cj^' 
FOR A OTIT-STEP ANGLE OF ATTACK 


(a) pi/c = 5/7 


y/i 





For n/l 

equal to - 





0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

^ig'fy/Wi) 

-1.3 

0.025339 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.2 

.067298 

.023339 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.1 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

0 

0 

0 

-1.0 

.198495 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

0 

0 

-9 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

0 

-.8 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

-.7 

.494931 

.622469 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

0 

-.6 

.494931 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

-.5 

.735859 

.622469 

.494931 

.382447 


.198493 

.126173 

.067298 

.023339 

0 

0 

-.4 

.930849 

.766859 

.622469 

.494931 

.382447 

.283 ©4 

.198493 

.126173 

.067298 

.023359 

0 

-3 

1.04930 

.897717 

.733726 

.589337 

.461799 

.349315 

.250723 

.165362 

.©3041 

.034167 

0 

-.2 

1.24247 

1.022886 

.834923 

.670937 

.526543 

.399010 

.28© 21 

.187932 

.102572 

.©5611 

0 

-.1 

1-42537 

1.157726 

.938137 

.750178 

.586189 

.441798 

.314260 

.201778 

.107498 

• 037253 

0 

0 

1.71432 

1.320358 

1.052711 

.833126 

. 6451 © 

.481175 

.3367© 

.213062 

.113216 

.©9149 

0 

.1 

1.98515 

1.589087 

1.195130 

.927484 

.707896 

.519936 

.359491 

. 22©04 

.119961 

.0413© 

0 

.2 

2.10430 

1.836762 

1.442695 

1.048740 

.781092 

.564900 

.587676 

.242903 

.128081 

. 0440 © 

0 

.3 

2.15443 

1.934844 

1.667303 

1.273240 

.882613 

.625377 

.424091 

.263557 

.1381© 

.151004 

.047229 

0 

.4 

2.14681 

1.958854 

1.739267 

1.474968 

1.091348 

.715331 

.474109 

.290757 

.© 12 © 

0 

.5 

2.08438 

1.920391 

1.735864 

1.526792 

1.277190 

.909457 

.550858 

.329142 

.. 16©92 

.©6609 

0 

.6 

1.96398 

1.823206 

1.670192 

1.500869 

1.308059 

1.076031 

. 7275 © 

.390548 

• 193©8 

.064041 

0 

.7 

1.77946 

L663767 

1.539261 

1.403506 

1.252644 

1 . 0797 © 

.869634 

.545674 

.©7054 

.075502 

0 

.8 

1.51367 

1.423440 

1.327050 

1.223027 

1.109179 

.981991 

.©5096 

.654029 

. 3637 © 

.096919 

0 

.9 

1.11123 

1.049985 

.984929 

.915250 

.839798 

• 756©5 

.©3525 

•554592 

.417548 

.181891 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


* 




PCm e '(y/j>’lA) 






- 1.3 

- 0.010999 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 1-2 

-.029761 

-.010999 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 1.1 

-.052085 

-.029761 

-.010999 

0 

0 

0 

0 

0 

0 

0 

0 

- 1.0 

-.076012 

-.052085 

-.029764 

-.010999 

0 

0 

0 

0 

0 

0 

0 

-.9 

-.100087 

-.076012 

-.052085 

-.029761 

-.010999 

0 

0 

0 

0 

0 

0 

-.8 

-.123018 

-.100087 

-.076012 

-.052085 

-.029761 

-.010999 

0 

0 

0 

0 

0 

-.7 

-143537 

-.123018 

-1000&7 

-.076012 

-.0520© 

-.029761 

-.010999 

0 

0 

0 

0 

-.6 

-.160289 

-.143537 

-.123018 

-.100087 

-.076012 

-.0520© 

-.029761 

-.010999 

0 

0 

0 

-.5 

-.171703 

-.160289 

-.143537 

-.123018 

-.100087 

-.076012 

-.0520© 

-.029761 

-.010999 

0 

0 

-.4 

-.175804 

-.171703 

-.160289 

-.143537 

-.123018 

-.100087 

-.076012 

-.© 20 © 

-.029761 

-.010999 

0 

-.3 

-.154246 

-.160190 

-.156087 

-.144674 

-.127922 

-.107401 

-.084472 

-.06© 97 

-.©6470 

-.014145 

0 

-.2 

-.107098 

-.127428 

-.133372 

-.129270 

-.117857 

-. 101104 

-i 080584 

-.©7654 

-•©3579 

-.012198 

0 

-.1 

-.031475 

-.074893 

-.095225 

-.101170 

-.097067 

-.085653 

-.068902 

-. 04 © 8 l 

-. 0275 © 

-.010100 

0 

0 

.109560 

.003450 

-.039971 

-.060299 

-.066244 

-.062141 

-.050727 

-.©5792 

-.020862 

-.007813 

0 

.1 

.251367 

.143258 

.039147 

-.004272 

-.024602 

-.030546 

-.028130 

-.021767 

-.013514 

-.005301 

0 

.2 

.329556 

.286136 

.180025 

.073915 

.050495 

.008549 

.002149 

. 0©845 

.005251 

. 0024 © 

0 

.3 

.382066 

.361737 

.318317 

.212207 

.104509 

.0554© 

.028715 

.012702 

.004245 

.000698 

0 

.4 

.415838 

.409894 

.389565 

.344554 

.233860 

.120461 

.067324 

.©5135 

.015493 

. 00442 6 

0 

.5 

.433166 

.437268 

.429691 

.404701 

.353949 

.238191 

.120487 

.063932 

.029415 

.008936 

0 

•6 

.434189 

.443882 

.443121 

.429608 

.399434 

.344338 

.225199 

• l©26l 

.047922 

.014697 

0 

.7 

.415320 

.426825 

.430171 

.423910 

.405836 

. 3721 © 

.314800 

.194884 

.076094 

.022748 

0 

.8 

.369481 

.380991 

.386357 

•384655 

.374563 

.354045 

.319554 

.263366 

.147246 

.©6210 

0 

• 9 

.281716 

.291131 

•296477 

.297212 

.292601 

.281579 

.262471 

.©2288 

.184276 

.082284 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABU! X.- PARTIAL SECTION LIFT AND MOMENT COEFFICIENTS BOjg' AND 1 


FOR A UNIT-STEP ARSLE OF A3TACK - Continued 


(b) fSl/c = 1 


y/i 





Far T]/l 

equal to - 





0 



0.3 





0.8 


^^jjl 

&Cj e '(y/l,T]/l) 

-0.9 

0.038956 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.8 

.113297 

.038956 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.7 

.214524 

.U 3297 

.038956 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

-341391 

.214524 

.113297 

.038956 

0 

0 

0 

0 

0 

0 

0 

-5 

.494934 

.343391 

.234524 

.113297 

.038956 

0 

0 

0 

0 

0 

0 

-A 

.678079 

.494934 

.341391 

.214524 

.113297 

.038956 

0 

0 

0 

0 

0 

-.3 

.896307 

.678079 

.494934 

.341391 

.214524 

.113297 

.038956 

0 

0 

0 

0 

-.2 

U59858 

.896307 

.678079 

.494934 

.343391 

.214524 

.113297 

.03895 6 

0 

0 

0 

-.1 

1.491333 

1.159858 

.896307 

.678079 

.494934 

.341391 

.214524 

.113297 

.038956 

0 

0 

0 

2.000000 

3,491351 

3,159860 

.896307 

.678076 

.494934 

.341389 

.234524 

.113298 

.098956 

0 

.1 

2.453263 

3,944634 

3-435965 

3,104472 

.840921 

.622691 

.439545 

.286005 

.159137 

.057904 

0 

.2 

2.678063 

2.346570 

1.857924 

1.329272 

•997784 

.734229 

.516001 

.332855 

.179313 

.061607 

0 

.3 

2. 794595 

2.531045 

2.199552 

1.690903 

1,182254 

.850761 

.587210 

.368980 

.193375 

.066116 

0 

.4 

2.825815 

2.607584 

2.344034 

2.012541 

1-503895 

.995243 

.665753 

.407062 

.211406 

.071799 

0 

-5 

2.778311 

2.595163 

2.376932 

2.113382 

3,761889 

3,273240 

.771201 

.460798 

.235749 

.079250 

0 

.6 

2.649749 

2.496209 

2.313064 

2.094830 

3,$51283 

1.506444 

1,018592 

.546767 

.271373 

.089656 

0 

-7 

2.428463 

2.301597 

2.148057 

3,964909 

1-753704 

3,511699 

3,217489 

.763944 

.331875 

.105701 

0 

.8 

2.085961 

3,984736 

1.657870 

3,712237 

3,552851 

3,374788 

3,169334 

.915640 

.509296 

.135688 

0 

-9 

1.544320 

3,469978 

3,378901 

3,281350 

i.175722 

1.059575 

.928935 

.776427 

.584570 

.254648 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


-0.9 

-0.017908 

0 

0 

0 

cr 

0 

0 

0 

0 

0 

0 

-.8 

-.047438 

-.017808 

0 

0 

0 

0 

0 

0 

0 

0 

0 

'•7 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

0 

0 

0 

-5 

-.143536 

-.314061 

-.080856 

-.047458 

-.017908 

0 

0 

0 

0 

0 

0 

-.4 

-.165594 

-.343536 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

0 

-3 

-.175684 

-.363594 

-.143536 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

-.2 

-.167131 

-.175684 

-.165594 

-.143536 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

-.1 

-.127232 

-.167131 

-.175684 

-.165594 

-.343536 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

-.127212 

-.167131 

-.175684 

-.165594 

-.143536 

-.114062 

-.080858 

-.047438 

-.017908 

0 

.1 

-.152677 

.025464 

-.101748 

-.341666 

-.150219 

-.140129 

-.118071 

-.088596 

-.055393 

-.021977 

0 

.2 

.235037 

.195119 

.067905 

-.059306 

-.099224 

-.107778 

-.097688 

-.075629 

-.04615 6 

-.017213 

0 

•3 

.292379 

.283825 

.245907 

.116694 

-.010518 

-.050436 

-.058989 

-.048899 

-.030353 

-.011857 

0 

.4 

.331928 

.342019 

•333467 

.293547 

.166336 

.039122 

-.000794 

-.012545 

-.011913 

-.005684 

0 

•5 

.355743 

.377801 

.387891 

.379538 

.339419 

.212207 

.081914 

.033149 

.010505 

.001664 

0 

.6 

.363565 

•393039 

.415098 

.425187 

.416633 

.373614 

.’237672 

.096958 

.039652 

.010875 

0 

• 7 

.353146 

.386348 

.415823 

.437881 

.444699 

.427101 

.373510 

.229183 

.082769 

.023446 

0 

.8 

.319064 

•352484 

.385686 

.411475 

.423574 

.418971 

.392499 

.333068 

.186742 

• 043635 

0 

• 9 

.247091 

.276623 

.505315 

.326265 

.358357 

.339959 

.328655 

.300000 

.244267 

.110947 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE X.- PARTIAL SECTION LIFT AND MDMHPT COEFFICIENTS pcj e ' AND BOmg' 
FOR A UNIT-STEP ANGLE OF A 3 TACK - Continued 

(c) pi/c - 4/3 






Far tj/ 

equal tc 

- 





y/i 












0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 


pcz 6 -’(y/i,n/i) 

-0.7 

0.021016 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

.113298 

.021016 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-5 

.253952 

.113298 

.021016 

0 

0 

0 

0 

0 

0 

0 

0 

-.4 

.440646 

.253952 

.113298 

.021016 

0 

0 

0 

0 

0 

0 

0 

-.3 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

0 

0 

0 

0 

0 

-.2 

.978403 

.678076 

.440646 

.253952 

.H3298 

.021016 

0 

0 

0 

0 

0 

-.1 

1.370753 

.978405 

.678076 

.440646 

•253952 

.113298 

.021016 

0 

0 

0 

0 

0 

2.000000 

1.370753 

.978405 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

0 

0 

.1 

2 . 62924 6 

2.000000 

1-370753 

.978405 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

0 

.2 

3.021605 

2.629246 

2.000000 

1.370753 

.976395 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

.3 

3.292095 

2.991777 

2.599419 

1.970174 

1.340925 

.948567 

.648249 

.410818 

.224124 

.083470 

0 

.4 

3.397278 

3.159848 

2.859530 

2.467172 

1.837926 

1.208678 

.816320 

.516002 

.278571 

.095730 

0 

■ 5 

3.379164 

3.192470 

2.955040 

2.654721 

2.262363 

1.633118 

1.003870 

.611522 

.314331 

.105669 

0 

.6 

3.242058 

3.101401 

2.914707 

2.677276 

2.576958 

1.984600 

1.355355 

.729023 

.361831 

.119542 

0 

•7 

2.970118 

2.877841 

2.737187 

2.550492 

2.313062 

2.012737 

1.623319 

I.OI8592 

.442500 

.140936 


.8 

2.507915 

2.507915 

2.394621 

2.253963 

2.067269 

1.833051 

1.558846 

I.220853 

.679061 

.180917 

0 

.9 

1.817493 

1.817493 

1.817493 

1.704192 

I.567623 

1.412759 

1.258580 

1.035232 

.779423 

.339531 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

EX^' (y/ 1 , 4 /*) 

- 0.7 

- 0.009946 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.5 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

0 

0 

-..4 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

0 

-.3 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

-.2 

-.175340 

-165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

-.1 

-.145363 

-.175340 

-.165594 

-.175340 

-.134344 

-.092126 

-. 047439 

-.009946 

0 

0 

0 

0 

.0 

0 

-.145363 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

.1 

.145363 

0 

-.145363 

-.175340 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

.2 

.175340 

.145363 

0 

-.145363 

-.175340 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

.3 

.179710 

.189456 

.159479 

.141157 

-.131247 

-.161224 

-.151478 

-.120229 

-.078010 

-. 033323 

0 

.4 

.202250 

.233500 

.243246 

.213269 

.067906 

-.077456 

-.107343 

-.097687 

-.066437 

-.026062 

0 

.5 

.225508 

.267723 

.298974 

.308722 

.278743 

.153380 

.011982 

-.041960 

-• 033713 

-.014650 

0 

.6 

.245058 

.289745 

.331962 

.363213 

.372959 

.342982 

.197619 

.050866 

.010187 

-.000590 

0 

•7 

.259451 

.296943 

.341631 

.383848 

.415098 

.424844 

.393465 

.237671 

.073396 

.018193 

0 

.8 

.275268 

.275268 

.322707 

.367394 

.409611 

.439329 

.437969 

.388647 

.218809 

.047776 

0 

• 9 

.250073 

.250073 

.250073 

.297512 

.340250 

.368951 

.377847 

.360794 

.304348 

.139585 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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For t 1/ Z equal to - 


0.1 0.2 


0.4 0.5 




0.113298 

0 

0 

0 

0 

.541391 

.113298 

0 

0 

0 

.678076 

.541391 

.113298 

0 

0 

1.159860 

.678076 

.341391 

.113298 

0 

2.000000 

1.159860 

.678076 

.341391 

.113298 

2. 840139 

2.000000 

1.159860 

.678076 

.341391 

3.321923 

2.840139 

2.000000 

1.159860 

.678076 

3.658608 

3.321923 

2.840139 

2.000000 

1.159860 

3.886701 

3.658608 

3.321923 

2.840139 

2.000000 

4.000000 

3.886701 

3.658608 

3.321923 

2. 840139 

3.387921 

3.837921 

3.724623 

3.496530 

3.159844 

3.503889 

3.503889 

3.503889 

3.390596 

3.162503 

2.991140 

2.991140 

2.991140 

2.991140 

2.877837 

2.199262 

2.199262 

2.199262 

2.199262 

2.199262 

G 

0- 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.113298 

0 

0 

0 

.341391 

.113298 

0 

0 

.678076 

.341391 

.113298 

0 

1.159860 

.678076 

.341391 

.113298 

1.837921 

.997782 

.515998 

.179313 

2.344029 

1.503889 

.663755 

.211409 

2.313059 

1.831270 

1.018592 

.271373 

1.857870 

1.552848 

1.169134 

.509296 

0 

0 

0 

0 


P^'Cy/W*) 


CO 

I 

• 

9 

0 

0 

0 

0 

-.114062 

-.047438 

0 

0 

0 

-.165592 

-.114062 

-.047438 

0 

0 

-.167131 

-.165592 

-.114062 

-.047438 

0 

0 

-.167131 

-.165595 

-.114062 

-.047438 

.167131 

0 

-.167131 

-.165592 

-.114062 

.165595 

.167131 

0 

-.167131 

-.165565 

.114062 

.165592 

.167131 

0 

-.167131 

.047438 

.114062 

.165595 

.167131 

0 

0 

.047438 

.114062 

.165592 

.167131 

.067907 

.067907 

.115345 

.181969 

.233502 

.166336 

.166336 

.166336 

.213774 

.280398 

.249504 

.249504 

.249504 

.249504 

.296942 

.271625 

.271625 

.271625 

.271625 

.271625 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

CD 

0 

0 

-.114062 

-.047438 

0 

-.165595 

-.114062 

-.047438 

-.167131 

-.165592 

-.114062 

.067907 

-.099223 

-.097688 

.335467 

.166336 

-.000793 

.415099 

.416632 

.237672 

.385685 

.423574 

.392499 

0 

0 

0 
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TAB1E I.- PARTIAL SECTICH LIFT AND MOMENT COEFFICIENTS PCj 0 ' AND pc^' 
FOR A UNIT-STEP ANOLE OF ATTACK - Concluded 


(e) pi/c - 4 


y h 


For T[] X equal to - 


0 

0.1 

a 


0.4 

0.5 



0.8 

— 


P<=x '(y/W*) 



CO 

8? 

3 

d 

0 

0 

0 

0 

0 

0 

0 

0 

0 

BE 

.678075 

CO 

a 

3 

0 

0 

0 

0 

0 

0 

0 

0 

la 

2.000000 

.678075 

CO 

a 

3 

0 

0 

0 

0 

0 

0 

0 


3.321923 

2.000000 

.678075 

.113298 

0 

0 

0 

0 

0 

0 

.2 

3.886711 

3.321923 

2.000000 

.678075 

.113298 

0 

0 

0 

0 

0 

•5 

4.000000 

3.886701 

3.321923 

2.000000 

.678076 

CD 

0 

0 

0 

0 

.4 

4.000000 

4.000000 

3.886711 

3.321923 

2.000000 

.678075 

.113298 

0 

0 

0 

.5 

4.000000 

4.000000 

4.000000 

3.886701 

3.321923 

2.000000 

.678076 

.113298 

0 

0 

.6 

4.000000 

4.000000 

4.000000 

4.000000 

3.886711 

3.321923 

2.000000 

.678075 

.113298 

0. 

.7 

4.000000 

4.000000 

4.000000 

4.000000 

4.000000 

3.886701 

3.321923 

2.000000 

.678076 

.113298 

.8 

3.837920 

3. 837920 

3.837920 

3-837920 

3.837920 

3.837900 

3.724633 

3.159044 

1.837926 

.515998 

•9 

2.991145 

2.991145 

2.991145 

2. 991^45 

2.991145 

2.991145 

2.991145 

2.877837 

2.313059 

1.018592 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


P'TBg'Cy/^n/i) 


- 0.2 

-0.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.1 

-.165594 

-.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.165594 

-.047438 

0 

0 

0 

0 

0 

0 

0 

.1 

.165594 

0 

.165594 

-.047438 

0 

0 

0 

0 

0 

0 

.2 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

0 

0 

0 

.3 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

0 

0 

.4 

0 

0 

.047438 

.165594 

0 

-.165594 

CO 

§ 

■3 

f 

0 

0 

0 

• 5 

0 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

.6 

0 

0 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

.7 

0 

0 

0 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

■ .8 

.679044 

.067904 

.067904 

.067904 

.067904 

.067904 

.115556 

.233502 

.067905 

-.097688 

.9 

.249504 

.259504 

.249504 

.249504 

.249504 

.249504 

.249504 

.296940 

.414640 

.237672 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE U.- SECTION LIFT AND MOMENT COEFFICIENTS 


P°lp AND pc^ 


FOR A UNIT RATE OF ROLL 


PCl P ” 1 + a + 1 + a PCi P 0 . 

p % * TTT pc ^. + I T7 P0m P o 


y/i 

Pc la (y/l) 

pc ZpQ (y/l) 


P=mp 0 (y/Z) 

(a) pz/e = 5/7 

0 

3. 714314 

0.306768 

0.109560 

-O.H8325 

.1 

3.589100 

.047812 

.145258 

-.159784 

.2 

3.442694 

-.184467 

.180026 

-.201629 

.3 

3.273239 

-387793 

.212207 

-.241915 

.4 

3.077662 

-.558969 

.240010 

-.278412 

• 5 

2.851240 

-.693521 

.261463 

-.308527 

.6 

2.586448 

-.784948 

.273900 

-.328680 

.7 

2.270441 

-.823060 

.273661 

-.333866 

.8 

1.878086 

-.789604 

.254594 

-.315697 

-9 

1.344769 

-.641916 

.202991 

-.255769 

1.0 

0 

0 

0 

0 

(b) pz/c = 1 

0 

4.000000 

0 

0 

0 

.1 

3.944613 

-.343532 

.025464 

-.025974 

.2 

3.837918 

-.631772 

.067905 

-.070621 

• 3 

3.690902 

-.873882 

.116694 

-. 123696 

.4 

3.503889 

-1.068886 

.166534 

-.179641 

• 5 

3.273240 

-1.212206 

.212207 

-.233427 

.6 

2 . 99 H 35 

-1.295676 

.249503 

-.279443 

• 7 

2.642985 

-1.305516 

.272287 

-.310407 

.8 

2.199256 

-1.216158 

.271624 

-.315084 

• 9 

1.583274 

-.966580 

.229183 

-.270436 

1.0 

0 

0 

0 

0 

(c) PI/ 0 = 4/3 

0 

4.000000 

0 

0 

0 

.1 

4.000000 

-.400000 

0 

0 

.2 

4.000000 

-.800000 

0 

0 

• 3 

3.970171 

-1.169877 

. 014115 

-.014256 


3.837919 

-1.433309 

.067906 • 

-.069943 

• 5 

3.633116 

-1.616485 

,133380 

-.140049 

.6 

3.355351 

-1.716783 

.197619 

-.211452 

• 7 

2.991142 

-1.719543 

.249502 

-.271957 

.8 

2.507913 

-1.593429 

.275267 

-.305547 

• 9 

1.817491 

-1.260633 

.250073 

-.282582 

1.0 

0 

0 

0 

0 
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TABLE II.- SECTION LIFT AND MOMENT COEFFICIENTS pc 2 ^ AND pc^ 


FOR A UNIT RATE OF ROLL - Concluded 


p %“m^a + — ■% 


p s = -— p x + — *%„ 


y h 

P c Za (y/z) 

P %o (y/i) 

pCn^Cy/ 1 ) 

P^Cy/z) 


(d) f!l/c = 2 


0 

4. 000000 

0 

0 

0 

.1 

4.000000 

-.400000 

0 

0 

.2 

4.000000 

-.800000 

0 

0 

• 3 

4.000000 

-1.200000 

0 

0 

.4 

4.000000 

-1.600000 

0 

0 

•5 

4.000000 

-2.000000 

0 

0 

.6 

3.857921 

- 2.254846 

.067906 

-.069264 

• 7 

3.503723 

-2.286270 

.166556 

-.172989 

.8 

2.991137 

- 2.145407 

.249504 

-.264474 

• 9 

2.199274 

-1.707722 

.271624 

-.293354 

1.0 

0 

0 

0 

0 



(e) pz/c = 

* 4 


0 

4.000000 

0 

0 

0 

.1 

4.000000 

-.400000 

0 

0 

.2 

4.000000 

-.800000 

0 

0 

.3 

4.000000 

-1.200000 

0 

0 

.4 

4.000000 

-1.600000 

0 

0 

• 5 

4.000000 

-2.000000 

0 

0 

.6 

4.000000 

- 2.400000 

0 

0 

• 7 

4.000000 

-2.800000 

0 

0 

.8 

3.837920 

-5.056584 

.067906 

-. 068585 

• 9 

2 . 99 H 35 

-2.567274 

.249505 

-.256987 

1.0 

0 

0 

0 

0 
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TABLE IU. - SECTIONAL LEFT AND MOMENT COEFFICIENTS gc^ AND gc^ 

FOR A UNIT AILERON DEFLECTION 

a = 0.l4 


V* = 

c a/ c 


y/i 

Pc^Cy/l) 

PCmgCy/i) 


(a) pz/c = 

= 5/7 

0 

-0.400000 

-0.160000 

.1 

-.646627 

-.251783 

.2 

-.760302 

-.301080 

.3 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-.'800000 

-.320000 

.6 

-.800000 

-.320000 

• 7 

-.800000 

-.320000 

.8 

-.742863 

-.292771 

• 9 

-.570248 

-.217640 

1.0 

0 

0 


(b) pz/c 

= 1 

0 

-0.400000 

-0.160000 

.1 

-.701015 

-.274665 

.2 

-.800000 

-.320000 

• 3 

-. 800000 

-.320000 

A 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.6 

-.800000 

-.320000 

• 7 

-.800000 

-.320000 

.8 

-.800000 

-.320000 

• 9 

-.654649 

-.253370 

1.0 

0 

0 


1.0 

0.2 


y/i 

pc^Cy/z) 

PCnigfr/*) 


(c) pz/c = 

= 4/3 

0 

-0.400000 

-0.160000 

.1 

-.749212 

-.295998 

.2 

-.800000 

-.320000 

.3 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-. 800000 

-.320000 

.6 

-. 800000 

-.320000 

• 7 

-. 800000 

-.320000 

.8 

-.800000 

-.320000 

• 9 

-.726623 

-.285314 

1.0 

0 

0 


(d) pz/c 

^ 2 

0 

-0.400000 

-0.160000 

.1 

-.800000 

-.320000 

.2 

-.800000 

-.320000 

.3 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.6 

-.800000 

-.320000 

.7 

-.800000 

-.320000 

.8 

-.800000 

-.320000 

• 9 

-.800000 

-.320000 

1.0 

0 

0 





TABIE IV. - INFLUENCE COEFFICIENT MATRIX 


USED IN SAMPLE CALCULATIONS 



0 

0 

0 

0 

0 


00000000 00 
0.055872 0, 244l88 0.449516 0.561924 0.625772 0.657864 0.676728 0.687456 0.695972 0.696760 

0.016848 0.075584 0.191504 0.415712 0.658496 0.761040 0.829440 0. 868248 0.891900 0. 909180 

0.005112 0.022520 0.058068 0.125552 0.24 q 8 o 4 0.476928 0.705708 0.855560 0.914 688 0.972596 

O.OOI656 0.Q07272 0.018900 0.040896 0.081596 0.155412 0.290576 0.556400 0.799256 0.991008 




r\ xvOi axo r\ /-\h r\ r\ r\on£nA r\ a)i )i rv£)i 

\j* wrr^vn- 


a aQp;Ao)i 

W* V/\J 


n i n oa^okk nacH^nkn 

vi xw _/wv vi v» / ry ^ ■ v» yyy v 1 v 



2 


Figure 1.- Configuration considered in tlie aeroelastic analysis. Positive 
directions of displacements and velocities are indicated t>y arrows . 
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(c) Wing with unit-step angle-of-attack distribution. 
Figure 2.- Configurations considered in the aerodynamic analysis 
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n 



0 


.2 .4 

Pressure 


.6 .8 

ra tio , & 


1.0 


(c) M = 1.338. 



0 .2 .4 .6 .8 1.0 

Pressure ratio , 

*0 


(d) M = I.667. 


Figure 3*- Continued. 
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O 12 3 

Moch number, M 


Figure 4.- Variation of pressure ratio at reversal with Mach number. 



0 12 3 

Mach number, M 

Figure 5 .- Variation with Mach number of the rate of roll for the rigid 

wing. 
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